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(54) Abstract Title 

Signal arrival direction estimation 

(57) In a DS-CDMA communication system for performing a transmission by a spreading process using a 
code having a period longer than a symbol length, received signals from an array antenna 1-0 to HM-1) are 
subjected to a despreading process by despreading units 3-0 to 3-(M-1). A cross-correlation function with time 
lags of the received signals is obtained by inverse modulating despread output signals by a known symbol in 
an inverse modulator 8. An arrival direction of a received signal is estimated by an arrival direction estimation 
unit 6 based on the cross-correlation function. A beam former 7 combines the despread output signals 
depending on the arrival direction and outputs an output signal to a channel receiver. The arrangement 
provides an accurate estimate of the arrival direction of a desired signal regardless of reception level and high 
interference. 
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At least one drawing originally filed was informal and the print reproduced here is taken from a later filed formal copy. 
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ARRIVAL DIRECTION ESTIMATION METHOD USING AH ARRAY 
ANTENNA AND PS -CDMA RECEIVER UNIT USING THE METHOD 

5 

The present invention relates to an arrival direction 
estimation method using an array antenna for estimating an 
10 arrival direction of a desired signal. The present 

invention further relates to a DS-CDMA receiver unit in -a 
Direct Sequence Code Division Multiple Access (DS-CDMA) 
system. 

In a base station of a DS-CDMA communication system, 
15 signal reception is performed by utilizing an array antenna. 
Further, an arrival direction of a desired signal is 
estimated from signals received from the array antenna. By 
transmitting and receiving in a beam directed in the 
estimated arrival direction, it is possible to reduce both 
20 interference and required transmission power due to an 

increase in antenna gain. Therefore, it is desirable that 
the arrival direction of the desired signal is estimated 
accurately. 

Fig. 7 of the dixiip srying drawings is a diagram shewing a reaeiwr unit cf a bsse 
25 station using an array antenna. In particular, Fig. 7 shows 
array antenna elements 71-0 through 71-(M-1) forming a 
receiving antenna, reception processors 72-0 to 72-(K-l), 
despreading units 73-0 to 73-(M-l), beam forming unit 74, a 
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1 



channel receiver 75, an arrival direction estimation unit 76 
and a beam former 77. 

Received signals r (0> to r CM " 1J from the array antenna 
elements 71-0 to 71-(M-1) are input to the reception 
5 processors 72-0 through 72-(K-l) corresponding to the 

different channels. In addition, the reception processor 
72-0 carries out a despreading process with respect to the 
received signals by the despreading units 73-0 to 73-(M-l), 
which are input with a spreading code Co- The despreading 

10 units 73-0 to 73-(M-l) provide despread output signals x (0) 
to x**" 1 * to the arrival direction estimation unit 7 6 and the 
beam former 77 of the beam forming unit 74. The arrival 
direction estimation unit 76 estimates an arrival direction 
of a desired signal based on a cross-correlation function 

15 among the received signals received by adjacent antenna 
elements. 

The beam former 77 multiplies a weighting coefficient 
obtained by estimating the arrival direction by the arrival 
direction estimation unit 76 to the despreading output 
20 signals x (0> through x (hhl) . Further, the beam former 77 then 
combines the products to provide an output signal to the 
channel receiver 75. Therefore, received data of the kth 
channel is re-generated by the reception processors 72-0 to 
72-(K-l) corresponding to channels 0 through K-l. 
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Fig. 8 cf the asrnparying drarirxp is a diagran showing a pronaslyKxiTsideaBa bean 
forming unit 74 including the arrival direction estination unit 76. 
In particular, Fig. 8 shows the beam forming unit 74 
including an arrival direction estimation unit 76 and a beam 
former 77 having multipliers 78-0 to 78-3 and an adder 79 . 

Received signals x (0) (n) to x {3> (n) from the antenna 
elements 71-0 to 71-3 (refer to Fig. 1) are input to the 
arrival direction estimation unit 76 and the beam former 77. 
The arrival direction estimation unit 7 6 outputs to the beam 
former 77 weighting coefficients aO to a3 obtained by the 
estimation of the arrival direction of the desired signal. 
The weighting coefficients aO to a3 are respectively input 
to the multipliers 78-0 to 78-3 of the beam former 77 to be 
multiplied with signals x (0) (n) to x <3) (n) and then coherently 
summed by the adder 79. An output signal y(n) of the adder 
79 is then input to the channel receiver 75 (refer to Fig. 
7) . 

ELg. 9 of the aascnpar^dng dravinge is a diagram stowing a poc^nojsly-acnsddened arrival 
direction estimation unit. In particular, Fig. 9 shows an 
arrival direction estimation unit 76 including correlation 
calculators 81-0 through 81-2, an adder 82, and an array 
weighting coefficient calculator 83. A correlation 
calculator 81 is shown in the lower part of Fig. 9, wherein 
each of the correlation calculators 81-0 to 81-2 include a 
multiplier 84 and an averaging filter 85. Further, received 



signals x (m> (n) and x Cml> (n) f rom nth and (m+l)th adjacent 
antenna elements are input to the multiplier 84. In this 
case, one of the received signals is a complex conjugate 
(indicated by *), and a complex conjugate product between 
5 the received signals from the adjacent antenna elements 

[cross-correlation function of the zero lag) is obtained, A 
temporal average is also performed by the averaging filter 
85 so that a cross-correlation function K M is output. 



10 at a time t satisfying (n-l)T < t < nT is described by the 
following formula (1), where 0 i tm) for 1=0 to N-l in 
formula (1) is described by formula (2), N denotes a number 
of users, A i denotes a received amplitude of an ith user 
signal, c A (t) denotes a spreading code of the ith user, T* 

15 denotes a relative delay of the ith user, bi denotes an ith 
user transmission symbol, d denotes a distance between the 
antenna elements, 6i denotes an arrival angle of the ith 
user signal, N (m) (t) denotes a noise signal, T denotes a 
symbol length, and a multipath is not taken into 

20 consideration for the sake of convenience. 



A received signal r 



im) ( t ) f rom a m th antenna element 



N-l 

r lm) (t)= Z AiCi(t-Ti)b(t-Ti)exp(-j0 
i=0 



25 



(1) 



=27nndsin6i 



(2) 
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A despreading output signal x (k , (lM {n) which is obtained 
by subjecting the received signal r fm) (t) from the mth 
antenna element to a despreading process by the spreading 
code c ()0 (t) of the kth user is described by the following 
5 formula (3), where w i3 (n)in the formula (3) is described by 
the formula (4) . 
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=A k b k (n)exp(-j0 t,B V 



N-l 



15 



(n-l)T 



i-0 
i*k 



20 



Ci (t-ti)c**(t-T k )bi (t-Ti)exp(-j0 tB ') Jdt 
+N (n) (n) 



N-l i 

=A k b k {n)exp(-j0 (n) )+ ZAj[ ( 



25 



i=0 J 
i*k 



(n-l)T 



Ci (t-Tj c k * (t-T k }bi (t-ti) dt] exp (- jc 
+N ,n) (n) 



30 



SA k b k (n)exp(-j0 w )f 



35 



N-l 

I [Ai w fi (n)exp (- j0i ,n) ) ] +N (n> <n) 



i=0 
i*k 



(3) 



40 



5 



15 



25 



.nT 
(n-1) 

(4) 



nT 

Wijtnja \ ci(t-ti)c k *(t-T k )bi(t-Xi)dt 

J (n-l)T 



In addition, the correlation calculator 81 of Fig. 9 
calculates a correlation R k ,Bl between the received symbol 
10 x ,m) (n) of the mth antenna element and the received symbol 
x tB+1) (n) of the (m+l)th antenna element based on the 
following formula (5), where brackets "< >" indicate 
averaging and the terms inside the brackets "< >" are shown 
in the formula (6) . 



R k (m, =<x k ,n) (n)x k (m+1, *(n)> (5) 



x k ,ffil (n)x k ,B+1, *(n) = |Aj 2 |b[ 2 exp[j (0 k , " , -a lt t - l, ) ] 

N-1 

20 + Z kl 2 |w lk (n) | 2 exp[j(el ) c (a,, 

i=0 
i^k 

-0 k ,B+1) )]+N tBl *(n)N (n, *n 



— (6) 

The correlation R k (ml is then averaged for each of the 
antenna elements and a cross-correlation function R k among 
the adjacent antenna elements is obtained based on the 
following formula (7) , where M denotes the number of antenna 
30 elements. 



M-2 

R k =[l/(M-l)]Z R k (n ' (7) 

m=0 



35 
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Fig. 9 shows a structure for adding the correlations 
R l0) to R {2> from the correlation calculators 81-0 to 81-2 by 
the adder 82. However, no calculation function related to 
the term [1/(M-1)] in the formula (7) is shown. Since (M-l) 
5 is a fixed value, it is clear that an averaging process can 
be carried out based on the addition. 

In formula (6) as described above, the second term on 
the right hand side of the equal sign is an uncorrelated 
signal and the term |w iIe (n) I is a small value. Hence, if 

10 there is no significant difference among the received signal 
amplitudes Ai of each of the users, it is possible to 
extract by averaging a phase error of a desired signal (kth 
user) of the first term on the right hand side of the equal 
sign of the formula (6) . Accordingly, based on a ratio of 

15 an imaginary part Im(R k ) and a real part Re (R k ) using the 
cross-correlation function R k of the formula (7), it is 
possible to estimate, from the following formula (8), an 
arrival angle 8 k of the kth user signal. 

20 e k =sin' 1 [{l/2ud)tan- 1 {Im(R k )/Re(R k ) }] (8) 

According to the previously-considered arrival direction 
estimation method described above, the arrival direction is 
estimated from the strongest signal received after the 

25 despreading process, regardless of whether or not the 

desired signal exists. For this reason, it is difficult to 
estimate the correct arrival direction of the desired signal 
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in cases where the reception amplitude of the desired signal 
is small and in cases where large interference signals 
exist. 

For example, when a multi-rate transmission including a 
5 high-speed rate and a low-speed rate in the DS-CDMA 

communication system is made, the received signal of the 
high-speed rate user becomes an interference signal for the 
received signal of the low-speed rate user since the 
transmission power of the high-speed rate signal is large. 
10 Consequently, a problem exists in that it is difficult to 
estimate the arrival direction of the desired signal of the 
low-speed rate user. 

15 

Accordingly, it is desirable ■ — 

to enable estimation of the correct arrival direction of the 
desired signal in a DS-CDMA communication system using a 
code having a period longer than the symbol length of the 
20 spreading code, even when a interference signal exists. 

An embodiment of e first aspect 

- of the present 

invention uses an array antenna in a DS-CDMA communication 
system for performing a transmission by performing a 
25 spreading process including a code having a period longer 
than a symbol length. The method also includes to subject 
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received signals from the array antenna to a despreading 
process. In this enbodine* a crtss^xxralatim functicn of lagp, preferably ecliriirg a zero 

time lag, for signals corresponding to adjacent antenna 

elements is obtained by inverse-modulating despread output 
5 signals by a known symbol. Further, an arrival direction of 
a desired signal based on the cross-correlation function is 
estimated. 

In another entedinent it is also possible to ise sn average value cf the 

cross-correlation function for a plurality of different lags 
10 as the cross-correlation function of the lags among the 
inverse modulated signals corresponding to the adjacent 
antenna elements, wherein the zero lag is excluded from this 
function. It is also possible to obtain the cross- 
correlation function using both positive and negative lags 
15 excluding the zero lag among the inverse modulated signals 
corresponding to the adjacent antenna elements. A pilot 
symbol may be used as the known symbol to inverse-modulate 
the received signals. It is also possible to provide a 
decision feedback symbol as the known symbol to inverse- 

20 modulate the received signals. 

An embodiment of a second aspect cf the present invention can provide 

a DS-CDMA receiver unit for use in a DS-CDMA cornnunication 
system, which carries out a transmission by performing a 
spreading process using a code having a period longer than a 
25 syTtal length. In this sitadiment the IB<jm receiver unit includes de^ceading units for 
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subjecting received signals from array antenna elements to a 
despreading process. Further included in this embodiment is a beam forming 

unit 4 for inputting despread output signals from the 
despreading units and a channel receiver for inputting 
5 output signals of the beam forming unit. The beam forming 
unit includes an inverse modulator for inverse-modulating 
the despread output signals by a known symbol, an arrival 
direction estimation unit for receiving signals from the 
inverse modulator and a beam former. The beam former 

10 provides an output by multiplying a weighting coefficient 
from the arrival direction estimation unit to the despread 
output signals and then adding such products. 

The inverse modulator may be configured to use either a 
pilot symbol or a decision feedback signal, as the known 

15 symbol. The arrival direction estimation unit may include a 
correlation calculator for obtaining a cross-correlation 
among output signals of the inverse modulator corresponding 
to the adjacent antenna elements and an array weighting 
coefficient calculator for obtaining and averaging the 

20 cross-correlation functions with respect to a plurality of 
lags excluding a zero lag. The array weighting coefficient 
calculator then outputs a weighting coefficient obtained as 
a result of estimating the arrival direction based thereon. 
In this case, the cross-correlation function with respect to 
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the plurality of lags may be calculated by providing first 

and second shift registers. 

Reference will now be made, by way of example only, to 

the accompanying drawings, in which: 

5 

Fig- 1 is a diagram showing a DS-CDMA 

receiver unit embodying the present invention; 

Fig. 2 is a diagram showing a spreading modulator of a 
transmission unit for use with an embodiment of the present invention; 
10 Fig. 3 is a diagram showing a despreading unit 

for use in an embodinent of the present invention? 

Fig. 4 is a diagram showing a beam forming unit 
for use in an entoodiment of the present invention; 

Fig. 5 is a diagram showing an arrival direction 
15 estimation unit for use in an embodiment of the present invention; 

Fig. 6 is a diagram showing a correlation calculator 
for use in an embodiment of the present invention; 

Fig. 7, discussed hereirbefeze, is a diagcan showing a receiver init of a boss 
station using an array antenna; 
20 Fig. 8, discussad bereirtefcre, is a diagram shewing a predajslynxreddened bean fanning 

unit having an arrival direction estimation unit; and 

Fig. 9, di an, mart hereirbefcre, is a diagran showing a p^vioBly-ocnsidQHi arrival 

direction estimation unit. 
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Fig. 1 is a diagram showing an embodiment of a DS-CDMA 
receiver unit according to the present invention. As can be 

5 seen from Fig. 1/ the receiver unit includes M array antenna 
elements 1-0 to 1- (M-l) forming a reception antenna, K 
reception processors 2-0 to 2-(K-l) provided in 
correspondence with K channels, M despreading units 3-0 to 
3- (M-l), a beam forming unit 4 and a channel receiver 5. 

10 The beam forming unit 4 includes an arrival direction 

estimation unit 6, a beam former 7, and an inverse modulator 
8. The reception processors' 2-0 to 2-(K-l) each have the 
same structure and receive signals from the array antenna 
elements 1-0 to l-(M-l). 

15 Fig. 1 further shows the receiver unit of the DS-CDMA 

communication system using a code C 0 having a period longer 
than a symbol length. Received signals r t0J to r (M_1) from the 
M array antenna elements 1-0 to l-(M-l) are respectively 
input. to the despreading units 3-0 to 3 (M-l) in the 

20 reception processor 2-0. For example, the despreading is 
carried out by the use of a spreading code Co. 

In this case, a transmission unit of a mobile station 
or the like would include a spreading modulator for the 
code. Fig. 2 shows an example of such a spreading 

25 modulator, which includes multipliers (exclusive-OR 
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circuits) 11, 12 and 14 to 17, a scramble modulator 13, 
polarity converters 18 to 21, and adders 22 and 23. 

With respect to a data symbol, a short code SCi of the 
quadrature code is multiplied by the multiplier 11 to 
provide code spreading. With respect to a pilot symbol/ a 
short code SC q of the quadrature code is multiplied by the 
multiplier 12 to provide code spreading. As a result, a 
distribution is made to quadrature channels (I,Q) . 

In addition, codes LCi and LC q of the quadrature code 
which have a period longer than the symbol length are also 
distributed to the quadrature channels (I,Q) and input to 
the multipliers 14 through 17. For example with respect to 
the data symbol which is subjected to code spreading by the 
short code SCi, the multiplier 14 multiplies the code LC 4 
and the multiplier 17 multiples the long code LC q to carry 
out a scramble modulation. With respect to the pilot symbol 
which is subjected to code spreading by the short code SC q , 
the multiplier 15 multiplies the code LCq and the multiplier 
16 multiplies the long LCi to also carry out scramble 
modulation. In this case, the pilot symbol may be inserted 
and transmitted for every predetermined length of data 
symbols, as a known symbol having a number of symbols less 
than a predetermined length. Alternatively, it is possible 
to transmit the pilot symbol together with the data symbol 
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having a phase that is orthogonal with respect to the data 
symbol. 

Output signals of each of the multipliers 14 to 17 are 
converted to the polarity converters 18 to 21. For example, 
5 a logic "0" is converted into "+1" or a logic n l" is 

converted into "-1"/ which are then input to the adders 22 
and 23. The adder 22 adds an output signal of the polarity 
converter 18 with a polarity and adds an output signal 
of the polarity converter 19 with a polarity. In other 

10 words, the adder 22 carries out a subtraction process. In 
addition, the adder 23 sums the output signals of the 
polarity converters 20 and 21. Accordingly, output signals 
of the adders 22 and 23 take a ternary value of +1, 0 or -1. 
The ternary value signals from the adders 22 and 23 are 

15 transmitted after carrying out a quadrature modulation or 
the like in a radio transmitter, which is not shown in Fig. 
2. The radio transmitter transmits signals subjected to 
spreading modulation by the code. 

In the embodiment of the DS-CDMA receiver unit shown in 

20 Fig. 1, the despreading units 3-0 to 3(M-1) carry out 
processes complementary to the process of the spreading 
modulator shown in Fig. 2. As shown in Fig. 3, for example, 
each despreading unit includes a descrambler demodulator 31 
and despreading units 32 to 35. The descramble demodulator 



31 includes multipliers 36-1 to 36-4 .adders 37-1 and 37-2. 

As can be seen from Fig. 3, all of the inputs and 
outputs include I and Q channel signals. However, in Fig. 
1, such detail was omitted. Instead, in Fig. 1, each of the 
5 I and Q channel signals are represented by a single signal. 
For example, receive : signal r(0), code C 0 , data signal D(0) 
and pilot signal x(0) of Fig. 1 correspond to I and Q 
channel signals, codes Ld and LCq, data signals Dj and D Q , 
and pilot signals x z and x Q , respectively in Fig. 3. 

10 The received I and Q channel signals are first 

subjected to the quadrature modulation in the radio 
receiver, the illustration of which is omitted. The same I 
and Q channel signals are then input to the descramble 
demodulator 31 to be descrambled using the codes LCi and 

15 LC q . The descrambled signals are then input to the 

de spreading units 32 to 35. Two of the descrambled signals 
are subjected to despreading using the short code SCi in the 
despreading units 32 and 33, which respectively output data 
symbols Di and Dq of the I and Q channels. While the other 

20 two descrambled signals are subjected to despreading vising 
the short code SC q in the despreading units 34 and 35, which 
respectively output pilot symbols x : and x Q of the I and Q 
channels . 
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As described above, the despreading units 3-0 to 3- (M- 
1) of Fig. 1 are configured to output the data symbols and 
the pilot symbols by carrying out the despreading process. 
In the present invention, the cross-correlation function 
which is used to estimate the arrival direction is obtained 
using the known symbol. This known symbol may be the pilot 
symbol, as described above. Alternatively, it is also 
possible to feed back a data symbol to the inverse modulator 
8 and to use this data symbol as the known symbol. In 
feeding back this data symbol, it would be judged by a 
judging unit included in the channel receiver. Fig. 1 shows 
the case where the pilot symbol is used as the known symbol. 

When signal x k tB,) (n) subjected to the despreading is 
further subjected to inverse modulation using a known 
transmission symbol bk(n) of the kth user, a symbol is 
denoted by z k (m) (n) , which is described by the following 
formula (9) . The formula (9) is developed into the formulas 
(10) and (11), where W i3 (n) in formula (11) is defined by 
the following formula (12) . 



z k w <n)«x k lm) (n)bk* (n) 



(9) 




(n-l)T i=0 
i*k 



N-l 

[ IAiCi(M0 



c k * (t-tk) exp (-j 0 (B ° ) ] dt+N w (n) 
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N-l JiT 

=A k exp(-j0 r "V IAj[ ( Ci(t-ti) 
i=0 J(n-1)T 

i*k 

Ck* (t-xOdtJexpC-j^V^n) 

(10) 



N-l 

z k ,^,) (n) 5 A k exp(-j0 (B, )+ I [A iWi3 (n) exp (-j0 i (B) ) 

i-0 
i*k 

+N (n "(n) (11) 

w i3 (n)s \ Ci (t-tjcic * (t-T k )dt 

J (n-l)T (12) 

In formulas (9) through (12), Ak denotes a reception 
amplitude of the kth user, c k (t) denotes a spreading code of 

the kth user, T k and Xi respectively denote relative delays 
of the kth and ith users, N'™ 1 (n) denotes a noise signal, and 
T denotes a symbol length. 

As may be seen from the following formula (13), a 
correlation R k tn) between a reception symbol z tBl (n) 
corresponding to the mth antenna element and L symbols 
before and after the reception symbol corresponding to the 
(m+l)th antenna element, is obtained by multiplying (1/2L) a 
summation of complex conjugate products amounting to 2L 



symbols in the positive and negative directions from -L to 
+L, excluding the zero lag, and averaging. The terms 
following the summation £ are shown in formula (14). 

L 

5 RK (m) =<U/2L) Z z [m> (n)2 ,n+1, *(n+l)> 

1=-L 

1* 0 (13) 

L 

10 Z z (B) (n)z ,m * 1> *(n+l) 

1=L 
1*0 

L N-l 

15 =|A k | 2 exp[j(0 lt ( » , -0 )t ,nfl) )]+ Z I|Ai| 2 W ik (n) 

1=-L i=0 
1*0 i*k 

20 W ik + (n+1) exp[j (0 i (B> -0 i ,m+1 >) ] 

L 

+ Z N <n " (n)N (n) *(n+l) (14) 

1— L 
25 1*0 

The formulas (13) and (14) show a case where both the 

positive and negative lags of the lags (time delays) other 

than the zero lag are used, where the correlation is 

30 obtained for the L symbols with respect to each of the 
positive and negative lags. When obtaining the cross- 
correlation for a plurality of different lags (time delays) , 
it is possible to obtain the correlation for the L symbols 
of a plurality of lags in only the positive direction or 

35 negative direction. In addition, it is also possible to 
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obtain the correlation for different numbers of symbols of 
different lags in the positive and negative directions. 

As described above, the correlation is obtained for L 
symbols of the different lags in each of the positive and 
5 negative directions, excluding the zero lag (time delay), of 
the reception symbol • z (ro) (n) corresponding to the mth antenna 
element and the reception symbol z <a * a) (n) corresponding to 
the (m+l)th antenna element. The arrival direction is 
estimated from this correlation. Hence, the arrival 
10 direction can be estimated by suppressing undesirable 

effects, even in cases where a phase deviation occurs due to 
fading. 

By obtaining the average of R <m) by formula (13), a 
cross-correlation function R* among the adjacent antenna 
15 elements can be obtained from the following formula (15) . 
In addition, an arrival angle 0* of the kth user can be 
obtained from the following formula (16) which is similar to 
formula (8) as described above, using the cross-correlation 
function R k . 

20 

M-2 

Rk-{1/(M-1) J X Rk (m) (15) 

m=0 

25 e^sin" 1 ! {l/2icd)tan" l {Im(Rk)/Re(R k } }] 

(16) 

Fig. 4 is a diagram showing the beam forming unit for 
this embodiment of the present invention. As can be seen 

19 



from Fig, 4, the beam forming unit 4 includes an arrival 
direction estimation unit 6, an inverse modulator 8 and a 
beam former 7. The beam former 7 includes multipliers 41 to 
44 and an adder 45. 

The despread output signals x t0) (n) to x (3) (n) 
corresponding to the A antenna elements are input to the 
beam forming unit 4, which are inverse modulated by the 
inverse modulator 8 using the known transmission symbol such 
as the pilot symbol, as described above in conjunction with 
formula (10) • As a result, inverse modulated signals zO to 
z3 are output from the inverse modulator 8. For example, 
the signals zO through z3 can be obtained by multiplying the 
pilot symbol by the despread output signals x {0) (n) to x (3) (n) 
upon reception of the pilot symbol, so as to obtain signals 
corresponding to the carriers of the desired signal 
excluding the pilot symbol component. 

The arrival direction estimation unit 6 calculates 
weighting coefficients aO to a3 based on the inverse 
modulated signals zO to z3, and inputs the weighting 
coefficients aO to a3.to the multipliers 41 to 44 of the 
beam former 7. As shown in Fig. 5, the arrival direction 
estimation unit 6 includes correlation calculators 51 to 53, 
an adder 54 and an array weighting coefficient calculator 
55. As shown in Fig^ 6, the correlation calculators 51 to 
53 are each provided with a first shift register 61, a 
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second shift register 62, an adder 63, a multiplier 64 and 
an averaging filter 65. 

The arrival direction estimation unit 6 shown in Fig. 5 

has a structure similar to that « - shown in 

5 Fig. 9. However, here, the correlation calculators are 

constructed differently. In the correlation calculator of 

this ennbodiment of the present invention, the inverse modulated output signal 

z (m) (n) corresponding to the mth antenna element is input to 
the first shift register 61, and the inverse modulated 

10 output signal z (m+X) (n) corresponding to the (m+l)th antenna 
element is input to the second shift register 62. The adder 
63 sums symbols L in both the positive and negative 
directions, excluding the symbol corresponding to the zero 
lag from the second shift register 62. An output of this 

15 adder 63 is then input to the multiplier 64. 

Accordingly, the multiplier 64 obtains a complex 
conjugate product of the inverse modulated output signal 
z (m> (n) corresponding to the mth antenna element obtained via 
the first shift register 61 and the inverse modulated output 

20 signal z (nr * 1} {n) corresponding to the (m+l)th antenna element 
for the L symbols in both the positive and negative 
directions excluding the zero lag at the center of the 
second shift register 62. An output signal of the 
multiplier 64 is then averaged by the averaging filter 65. 

25 In other words, it is possible to obtain the correlation R (n0 
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from the following formula (17) . In addition, formula (18) 
is obtained by developing the terms following the summation. 
Further, by averaging the correlation R lBl of the formula 
(17) for each of the antenna elements, the cross-correlation 
function R* among the adjacent antenna elements can be 
obtained from formula (19) . 



L 

R k ,B, =< Z z (B " (n)z <m+1) *(n+l)> 
1— L 

1* 0 (17) 



L 

Z 2 <B " (n) z 1 -"* (n+1) 
1=-L 

1*0 



=|A k | 2 exp t j (0 k ,B) -0 te (0+1 ' ) ] +Z 2>i| 2 W ilt (n) 

1=-L i=0 

1*0 i*k 

W ik * (n+ 1 ) exp [ j ( 0 4 M -0i tm * 1> )) 
L 

+ Z N (,B) (n)N ,,B, *(n+l) — - (18) 

1=-L 
1*0 

M-2 

R k ={l/(M-1) } Z R k tB " (19) 

m=0 

In Fig. 5, the output signals R 110 ' through R <2) of the 
correlation calculators 41 to 53 are then summed by the 
adder 54. However, the term [1/(M-1)] of formula (19) 
averages depending on the number of array antenna elements, 
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which is determined in advance- Accordingly, the function 
of carrying out the process that corresponds to the term 
[1/(M-1)] is not included in Fig- 5. 

Using the cross-correlation function R k of the formula 
(19), the arrival angle 6 k of the kth user is obtained from 
the following formula (20) . 

9 k -sin" 1 ! (l/2icd)tan- l {Im(Rk) /Re (R k ) }] 

(20) 

That is, the arrival angle 6 k is obtained based on the 
ratio of the imaginary part Im(R k ) and the real part Re (R k ) 
of the cross-correlation function R k . 

The weighting coefficient a k (m} input to the beam former 
7 can be obtained from the following formula (21) . Further, 
output signals y k ,o(n) and y k , P (n) of the beam former 7 are 
described by the following set of formulas (22) ♦ 

a k (m> =(R k *) ra (m=0 to M-l) (21) 

M-l 

yk,o(n)= I a k Cn> D k {m) (n) Data Signal 

m=0 

M-l 

y*, P (n)= Z a k {a, x k <m) (n) Pilot Signal 

m=0 

(22) 

In other words, in Fig. 4 the weighting coefficients 
a (0> to a <3) from the arrival direction estimation unit 6 are 
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obtained from formula (21) . The weighting coefficients a (0> 
to a (3) are respectively input to the multipliers 41 to 44 of 
the beam former 7 and multiplied by the despread output 
signals x (0> (n) through x (3> (n) . The adder 45 then produces an 
5 output y{n) by summing the products of the inverse despread 
output signals x (0} (n) to x (3} (n) and weighting coefficients 
a (0> to a (3) while matching the phase with the estimated 
arrival direction. The phase is matched with the arrival 
direction because the weighting coefficients a (0) through a <3> 

10 correspond to the complex conjugates of the cross- 
correlation function R k . The output y(n) of the adder 45 is 
then input to the channel receiver 5 (refer to Fig, 1), so 
that the pilot signal y P and data, signal y D corresponding to 
the channel are demodulated and output. 

15 It should be noted that the present invention is not 

limited to the embodiment described above, and various 
variations and modifications may be made. For example, the 
functions of the operation processes of the arrival 
direction estimation unit 6 and the like can be realized by 

20 the functions' of a digital signal processor (DSP) , and other 
functions can similarly be realized by the functions of the 
processor. In addition, it is possible to provide an 
interference canceller or the like to eliminate the 
interference among the channels. 



As described abcxe, aooauling to an ertodinEnt of the present inv^rrticri, 
a DS-CDMA communication system performs a transmission by a 
spreading process using a code having a period longer than a 
symbol length. Further , received signals from corresponding 
array antenna elements 1-0 to l-(M-l) are subjected to a 
despreading process. A cross-correlation function is 

also obtained for lags (preferably excluding the zero lag) for signals 
corresponding to adjacent antenna elements obtained by 
inverse-modulating despread output signals by a known 
symbol. An arrival direction of a desired signal is 
estimated based on a cross-correlation function. This 
erTtoodiinent has an advantage in that the arrival direction 
of the desired signal can be estimated accurately, even in 
cases where the reception level of the desired signal is low 
or the level of the interference signal is high. Therefore, 
it is possible to increase the antenna gain of the array 
antenna, to eliminate the interference, and to reduce the 
required transmission power. 

In another enbodiment, by obtaining the cross-correlation function 

for a plurality of symbols in the positive and negative 
directions other than the zero lag among inverse modulated 
signals corresponding to the adjacent antenna elements and 
averaging the cross-correlation functions, there is an 
advantage in that the degree of estimating the arrival 
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direction of the desired signal can be improved, even with 
respect to a phase deviation caused by fading or the like. 
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CLAIMS ; 

1. An arrival direction estimation netted for use in a DB-Cim 
communication system including an array antenna for 
performing a transmission by a spreading process using a 

5 code having a period longer than a symbol length, said 
method comprising the steps of: 

subjecting received signals from the array antenna 
to a despreading process in order to produce despread output 
signals; 

10 obtaining a cross-correlation function of time 

lags of the received signals from adjacent antenna elements 
by inverse-modulating the despread output signals by known 
symbols; and 

✓ 

estimating an arrival direction of one of the 
15 received signals based on the cross-correlation function. 

2. The method of claim 1, wherein said cross- 
correlation function of the time lags between 'the received 
signals of adjacent antenna elements exclude a zero time 
lag. 

20 3. The method of claim 1 cr 2, further comprising the steps 

of: 

averaging values of the cross-correlation function 
for a plurality of time lags as the cross-correlation 
function of the time lags of the received signals from the 
25 adjacent antenna elements. 
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4 . The arrival direction estimation method of any preceding claim f 
wherein said obtaining the cross-correlation function uses 
both positive and negative time lags between received 
signals. 

5 5. The method of any preceding claim, wherein the known symbol is 

a pilot symbol. 

6 . The method of any one of claims 1 to 4, further ccnprising the step 

of: 

providing a decision feedback symbol as the known 

10 symbol. 

7. An arrival direction estimation method for use in a DS- 
CDMA communi cation system using signals subjected to a 
spreading process by a code, said method comprising the 
steps of: 

15 subjecting the signals to a despreading process in 

order to produce despread signals; 

inverse modulating the despread signals by a known 
symbol in order to obtain inverse modulated signals; and 

producing weighting coefficients corresponding to 
20 arrival directions of the signals from the inverse modulated 
signals. 

8. The method of claim 7, wherein the code has a period 
longer than a symbol length. 

9. The method of claim 7 or 8, wherein the producing of 
25 weighting coefficients includes obtaining a cross- 

28 



correlation function of time lags from the inverse modulated 
signals and producing the weighting coefficients based on 
the cross-correlation function. 

10. The method of any one of claims 7 to 9, wherein the cross- 
5 correlation of the time delays excludes a zero time lag. 

11. The method of any one of claims 7 to 10, wherein the knowi syrbol is 

a pilot symbol. 

12 . The nett ed of eny cneaf plains 7 to 10, wherein tte known synbaL is 
a decision feedback symbol obtained from a channel receiver. 
10 13. The method of any one of claims 7 to 12, which further includes 

multiplying the despread signals by the weighting 
coefficients to produce weighted signals. 

14. The method of claim 13, which further includes 
summing the weighted signals in order to produce output 

15 signals. 

15. A DS-CDMA receiver unit for use in a DS-CDMA ccmnunication 
system to produc e a transmission by a spreading process 
using a code having a period longer than a symbol length, 
said receiver unit comprising: 

20 a despreading unit for subjecting received signals 

to a despreading process to produce despread output signals; 

a beam forming unit for receiving said despread output 
signals from said despreading unit; and 

a channel receiver for receiving output signals 
25 from said beam forming unit; 
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said beam forming unit including an inverse 
modulator for inverse-modulating the despread output signals 
by a known symbol, an arrival direction estimation unit for 
producing weighting coefficients based on signals from said 
5 inverse modulator, and a beam former for producing an output 
by multiplying the weighting coefficients from said arrival 
direction estimation unit with the despread output signals 
producing weighted signals and then combining the weighted 
signals. 

10 16,. The IDS-CDMA receiver unit of claim 15, wherein said 

inverse modulator inverse-modulates the despread output 
signals from said despread unit using a pilot symbol. 

17. The D5-CDMA receiver unit of claim 15, wherein said 
inverse modulator inverse-modulates the despread output 

15 signals from said despreading unit by a decision feedback 
symbol from said channel receiver, 

18. The DS-CDMA. receiver unit of any one of claims 15 to 17, vtereinsarid 

arrival direction estimation unit includes a correlation 
calculator for obtaining a cross-correlation of the signals 
20 from said inverse modulator, and an array weighting 
coefficient calculator for calculating the weighting 
coefficients based on the cross-correlation obtained by said 
correlation calculator. 



19. The DS-CDMA receiver unit of claim 18, wherein said 
cross-correlation includes a plurality of positive and 
negative time lags excluding a zero time lag. 

20. The DS-CDMA receiver unit of claim 18 or 19, wherein said 
5 correlation calculator includes a first shift register and a 

second shift register for receiving signals from said 
inverse modulator, a multiplier for obtaining the 
correlation, and an averaging filter, said second shift 
register provides to said multiplier a plurality of positive 
10 and negative symbols with respect to a symbol input to said 
multiplier from said first shift register. 

21. An arrival direction estimation method substantially as 
hereinbefore described with reference to Figures 1 to 6 of 

the accompanying drawings. 

22. A CDMA receiver unit substantially as hereinbefore 

" described with reference to Figures 1 to ^ of the accompanying 
drawings. 
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